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ABSTRACT 


The NAVSTAR Global Positioning System is a worldwide, 
all-weather, satellite positioning system capable of high 
accuracy real-time poSition determinations. The Applied 
Research Laboratories (ARL:UT), in conjunction with the 
Naval Surface Weapons Center (NSWC/DL) andi other government 
agencies, conducted geodetic fieid tests of a jovernment 
sponsored prototype receiver, the TI4¥100/GEOSTAR, in March 
1984. Data were acquired from four satellites by three 
receivers, with antennas located at known Stations, over 
approximately six-hour periods each day. 

Based on the real-time solutions acquired on March 1, 8, 
and 9, 1984, absolute point position determinations have an 
average discrepancy of 7.5 meters with a one Sijma repeat- 
ability of less than a meter. Relative positions were 
determined to an average relative accuracy of 1:9,700 for 


distances of 14 to 2€ kilometers. 
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I. INTRODUCTION 

The Navai Surface Weapons Center (NSWC), with sponsor- 
Ship by the Defense Nappitg Agency (DMA), has developed a 
system for geodetic positioning with the use of the NAVSTAR 
Global Positioning System (GPS) [kef. 1]. A government 
sponsored prototype receiver, the TI4100/GSOSTAR, was devel- 
oped by Texas Instruments through contract with Applied 
Researcn Laboratories, The University of Texas at Austin 
(ARL:0T) [{Ref. 2]. The receiver was tested by ARL:UT in 
conjunction with NSWC and other government agencies over 
short and medium base lines in March 1984. The Defense 
Mapping Agency Hydrographic/Topographis Center (DMAHTC) nas 
expressed interest in independent analyses of the data from 
these tests [ Ref. 3}. 

Absolute point positioning, as it is used in this 
thesis, is the determination of a position with the use of 
one GPS receiver. Relative positioning is the determination 
of two positions relative to each other with the use of two 
receivers. The basic premise orf real-time relative posi- 
tioning is that GPS observation errors are to a large extent 
systematic in the same yeneral area at the same time. If 
such is the case, relative positioning by two receivers can 
locate positions relative to each other more accurately than 
absolute positioning alone. 

Accuracy standards for horizontal control require 
certain relative accuracies between directly cornected adja- 
GCeENe pormes. “Third-order surveys are used for local hori- 
zontal control in hydrographic surveying and other local 
projects. Third-order, Class II standards reguire relative 
accuracies of at least 1:5,000 between directly connected 
adjacent points (e.g., 5 meters relative accuracy between 


Stations separated by 25 kilometers). [{Ref. 4] 
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Geodetic positioning with GPS has advantayes over 
conventional surveying methods. Visibility between sites 
and fair weather conditions are not required. Real-time 
positioniny has a Special application in the case where the 
position of a station is needed in a Short amount of time 
and before the data can be post-processed. An analysis of 
the real-time solutions of the GEOSTAR short and medium base 
line tests will give an estimate as to what degree of accu- 


racy the system can achieve in real-time. 
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Ii. BACKGROUND 


Se ee eS ce 


eee LueeGLObAL POST? PONENG SYSTEM 


The Global Positioning System 1s a worldwide, all- 
weather, satellite positioning system. GPS can be thought 
of aS consisting of three components; the Space Segment, the 
Control Segment, and the User Segment. The Space Segment 
will consist of eighteen satellites (with three spares) in 
Six Orbit planes when fully operational in late 1988 or 
early 1989 [Ref. 5]. The satellites will orbit the Earth in 
twelve-hour periods at altitudes of approximately 20,000 
Kilometers. This constellation will provide four satellites 
in view at one time at elevations of greater than 20° above 
the horizon. {Ref. 6} The satellites" ground paths will be 
fixed. The time at which a particular satellite will be in 
view at a certain position will be about frour minutes 
earlier each day due to the differences between solar and 
Sidereal time. [Ref. 7] 

The GPS satelkites have on orbit weight of 950 vounds. 
Power is supplied by solar panels and nickel-cadmiun 
batteries. The satellites are equipped with atomic clocks 
Vitec litheseOnmenesoLaer-of I part in 1043. [Ref. 8] 
Perturbations in the satellite orbits are caused by the 
forces of attraction of astronomical bodies, gravitational 
irregularities or the earth, and solar radiation pressure. 
These perturbations are compensatei for by the Control 
Segment and frequency standard offsets. [Ref. 7] 

The Control Segment consists of four monitor stations 
locatea at Alaska, Hawaii, Guam, and Vandenburg AFB; anda 
Master Control Station (MCS) and an Upload Station at 
Vandenburg AFB {[Ref. 9J. The monitor stations track the 
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satellites then send the data tovuthe MCS tereproce-- 17 ee 
determine updated predictions of satellite and clock 
behavior. The updated predictions are transmitted to the 


Satellites by the Upload Station. [Ref. 10] 


B. TRANSMISSION SIGNALS 


The GPS NAVSTAR satellites simultaneously transmit navi- 
gation information on 2 RF frequencies; the L1 frequency at 
1575.8 MHz.§and the L2 at 1227.6.8HzZ. sobh TEsaiencie sea 
modulated with pseudo-random noise (PRN) codes. The Precise 
Positioning Service (PPS), formerly called the "P-code", is 
modulated on both freguencies and the Standard Positioning 
Service (SPS), formerly called the "C/A-code", is avaiiabie 
only on the L1 freguency. PRN codes, different for each 
Satellite, serve the purpose of uniquely identifying the 
particular satellite and provide a means of determining the 
Signal transit time by measuring the phase shift necessary 
zor the receiver to match codes. [Ref. 11 and 12] 

The SPS is a relatively short code of 1023 bits. The 
code repeats itself every millisecond at a bat rate can] 
Mbps. Due to the short duration time of the SPS, the code 
is relatively easy to match and lock on to. The PPS has a 
7-day code duration and 10 times the bit rate, making the 
PPS code relatively difficult to match with the receiver's 
internally generated code. The PPS offers the advantage of 
a finer resolution of the phase shift necessary to match 
codes, thus a greater measurement accuracy of signal transit 
tine. (PRetite ia 


C. THE NAVIGATION MESSAGE 


The GPS Navigation Message is modulated on both the L1 


and L2 frequencies and contains information necessary to 


compute the user's position. The information is in the fora 


12 


Ob ya 20m ePsemo meade sGercal COMMON tO DOLLA the PPS. and SPS. 

The data frame is 1500 bits iong and consists of 5 subfranes 
of 6 seconds (300 bits) each. The Navigation Message 
provides information such as satellite e,hemerides, system 
time, sateliite clock behavior, and transmitter status. 
PRefs 13 | 

Subframes are divided into ten 30-bit words. The first 
two words of each subframe contain system time and time 
synchronization information for transfer from the S2S to the 
PPS in the form of hand-over words (HOW) and telemetry words 
(TLM). The remaining 8 words contain user navigation data 
generated by the Control Segment. [Ref. 13] 

Data Block 1 occupies words 2 through 10 in subframe 
one. It contains frequency standard corrections, an associ- 
ated age of data word, and ionospheric delay corrections for 
Single frequency receivers. The clock correction parameters 
are used to determine the satellite's clock offset from GPS 


time by the following equations: 


t=t - At oe 1 


SV SV 


NE =" 205 a, (c— too) ios Ca Las (25.2) 
where t is the GPS time in seconds measured froin the begin- 
hing of the GPS week (approximately Saturday/Sunday midnight 
GMI) , ae is the satellite's PRN code phase time at tne tine 
of message transmission, t,. 1s Data Block 1 reference tine 
(seconds), and Anse Ais and a, are Data POG wesc LOCK. COLL oc— 
tion parameters (coefficients of a second order polynomial). 


y ih equation 2.2. 


Note that t may be approximated by t . 
fRef£. 13] 

The age of data word (AODC) is the time between Data 
Block 1 reference time and the tine since the last measure- 


went used in computing the clock correction parameters 
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[fRef. 13]. Test data in this thesis wete opraine si 
both the L1 and L2 frequencies, therefore, the 1:onospheric 
delay parameters for single frequency receivers contained ina 
Data Block 1 are not used (see Van Dierendonck, et al, 

137 oe 

Data Block 2 contains ephemeris information in subfranes 
2 and 3. This data includes associated age of data words 
(AODE) which represent the time difference between Data 
Plock 2 reference time and the time of the last measurement 
used in the estimation of the ephemeris parameters. fhe 
ephemeris data includes conventional Keplerian-type parame- 
ters, secular drift terms, and harmonic coefficients used to 
cocrect for perturbations in the Keplerian vorbit. | \iimet iS 
The satellite's position is calculated Ly Van Dierendonck's 
algorithm [{Ref. 13] with the addition of a secular drift 
term to correct for perturbations in the inclination of the 
satellite's orbit [{Ref. 14]. 

The Message Block appears in the third through the tenth 
30-bit words of subframe 4 of the Navigation Message. This 
space is reserved for alphanumeric messages provided to the 
user by the Control Segment. Space exists for 23 eight-bit 
ASCI#= characters. “Reb. 159) 

Almanac information appears in Data Block 3 in the fifth 
subframe. Data Block 3 consists of 25 subframes allowing 
for almanacs of up to 24 satellites (with one dumny 
subframe). Twenty-five navigation messages must be received 
to acquire the complete data set of Data Block 3. Extra 
subframes can be used to repeat the almanacs for certain 
Satellites. Satellite almanacs provide the user with less 
precise position and clock correction information. Almanacs 
aid the user in satellite selection and simplify the acoui- 
Sition of Signals. Almanacs include ephemeris parameters, 
Clock correction parameters, satellite identification, and 
satellite health terms. [Ref. 13] 
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Do SaeGEIVER OPERATION 


The GEOSTAR GPS receiver iS a governnent sponsored 
version of the TI4100 commercial receiver (Fig. 2.1). #£=The 
GEOSTAR 1S a sinyle channel, multiplex receiver capable of 
Peco wri Lolresateltates On ehe DPiwand LZ PPS. The control 
and display unit (CDU) is a hand-held input/output device 
With two l16-character alphanumeric displays. The antenna 
assembly consists of Guconi Cammonniagibectlohnal antenna with 
a preamp that allows remote operation of the receiver. Two 
100-foot cables connect the antenna assembly to the 
receiver. One cabdle is used to input DC power to the preamp 
and output the amplified RF Signals. The other cable is 
used for a built-in test signal injection from the receiver. 
{Ref. 15] 

A typical GEOSTAR geodetic surveying unit consists of a 
2-man team, a Station wagon or van, and the sSSOSTAR field 
set. The field set contains the receiver, CDJ, the antenna 
assembly, a cassette tape recorder, a tripod with a leveiing 
bracket, two 12-volt automobile batteries, an EFRATOM FRY-H 
rubidium oscillator, a weather station with sensors, cables, 
transportation cases, tools, manuals, and spare parts. The 
tield set 1s easily contained and transported in the rear of 
the vehicle. 

The antenna assembly iS mounted on a standard surveyor's 
tripod with a leveling bracket then leveled and plumbed over 
the mark to be surveyed. A reference mark on the base of 
the antenna assembly is oriented to north usSiny a hand 
compass. The height of the L2 phase center mark at the base 
of the antenna assembly above the survey point is measured 
and recorded. The antenna assembly, the power supply, tne 
cassette recorder, the weather station, and the rubidiun 
oscillator are then connected to the receiver with the 
appropriate cables. Scratch tapes are inserted in the 


cassette recorder for the power-up procedure. 
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The TIG100/SEOSTAR Receiver 


Figure 2.1 


The system is turned on and the receiver and cassette 
RecoOrpG@e: Seic-testS abe initiated. When the selr-tests are 
successfully completed, the GESAR bootstrap tape is inserted 
i Wetve ieee Drogmay 2S loaded ani blank tapes are 
iPNSeo eet etn eomewOomtape Crives. ~The operator is then 
Peete o so yetnes COU rom 1nput information tncluling an esti- 
Paco wOmmenemUSor Ss sDosition, weather data, Satellites to ke 
Pedekea, and Doppler alding values for each satellite. fThe 
system tnen operates unattended except ror the replacement 
of cassette tapes as required and the periodic manual input 
of weather data (1f the automatic weather station 1s not 


used). 


E. POSITION COMPUTATION 


The satellite's position and the time of signal trans- 
mission are computed with the broadcast epheaeris and the 
clock correction parazeters included in the Navigation 
Message. The time between the transmission and the recer- 
Prommonmere Sitetiite Sigmal multiplied by the speed of 
iigat is the pseudorange. 

P seudoranges used to calculate the real-time solutions 
are corrected for atmospheric delays in the GEOSTAR Receiver 
[Rete sli j= LomospheriG retraction iS nearly inversely 
proportional to the square of the freguency. By receiving 
coth the Li and the L2 frequencies, the ionospheric delay is 
calculated by comparison. [Ref. 16] Tropospheric delay is 
calculated by a model based on temperature, pressure, 
humidity, and the elevation of the satellite above the 
horizon [Ref. 17]. Knowing the positions of four satellites 
and the ranges to the user, the user's position is calcu- 


lated by a least Squares method. 
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F. GEODETIC DATUMS 


The GEOSTAR real-time solutions are referenced in DoD Worid 
Geodetic System 1972 (WGS-72) Cartesian coordinates. WGS-/2 
is an earth-centered, earth-fixed, glonpal, geodetic datum 
adopted by DMA in March 1974 [Ref. 18]. The reference 
surface is an ellipsoid that kest approximates the shape or 
the earth's geoid. The ellipsoid is normally defined by the 
length of the semimajor axis and the flattening. 

Latitude ranges between -$0° and 90° with north fosi- 
tive. Longitude is measured positive east from 09 to just 
less than 3609 from the WGS-72 reference meridian at 0".54 
east of the Greenwich Meridian. The Cartesian X-axis 1s 
formed by the intersection of the reference meridian plane 
and the plane of the equator. The Y-axis is 90° east of the 
X-axis in the plane of the equator. The Z-axis 1S in tne 
direction of the earth's mean rotation axis based on the 
average terrestrial pole of 1900-05. [Ref. 19] 

The North American Datum 1927 (NAD-27) is a regional 
datum based on the Clarke 1866 Hilipsoid with the origin at 
Meades' Ranch, Kansas. The Clarke 1366 Ellipsoid minimizes 
the deflection of the vertical and geoidal undulations in 
the continental United States. The National Ocean Survey 
(NOS), the National Geodetic Survey (NGS), and the United 
States Geological Survey (USGS) use NAD-27 coordinates in 
most survey operations. The transformation from WAD-27 to 
wGS-72 reference systems creates offsets of the order of 5 
meters [Ref. 2]. It 1S assumed that tails error is removed 
when the reverse transformation is made; but the error 
source will still exist in the analysis of the GS-72 point 


position discrepancies. 
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Tit. AKALYSIS SETHODS 


A. EXPERIMENT DESIGN 


Relative position field testing of the GEOSTAR Receiver 
waS performea in central Texas in February and March 1984. 
Three receivers, each with an antenna located over a 
different Known position, were operated continuously over 
approximately the same period of satellite visibility. Four 
satellites, with PRN codes 6, 8, 9, and 11, were tracked in 
these tests. The period of data gathering was coordinated 
with the time period in which all four satellites were in 
view at the same time. The real-time positions of the 
antennas were recorded every 12 seconds over approximately a 
Six-hour period each day. The data were recorded on 


cassette tape then transferred to 9-track tape by ARL: UT. 


B. GPS USER TAPES 


The recorded data and the associated position coordi- 
hates used in the analysis were obtained from NSYC/DL. The 
GPS User Tapes are a standardized format, 9-track magnetic 
tape. The tapes are unlabelled, ASCII coded, 1600 bDpi, with 
80 characters per physical record and 1 record per block. 

The first file on the tapes consists of header comment 
records. This file assists the user in identifying the tape 
and includes information such as the date and site of the 
data acqums itaon. 

Data files begin with an 80-character control record 
(Cata Item List) which identifies the types of data that 
follow. The first location of all records is reserved for a 
ComerpoOlecitede ten. ine Data Litem List control record has an 


aSterisx ‘*' in the first position followed by an integer in 
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the next three positions which specifies the numpoer of times 
to repeat the sequence of the Data Item List. The remaininy 
76 characters contain up tO” 19ueedeegit "data Cypeuticr en. 
fiers which specify the type Gewdata that 15) come oO eae 
control record. Control records ame redcdeneecus 
(AT 73 5104 jp ehOrtiac. 

Data information are one or more records in length and 
divided into groupings called @data tenis (lances, 
{[Ref. 19]. 


PARLE. © 
Data Types of GPS User Tapes 


item Identifier Data Type 

VOC T=soe8 Exchanje conversion programa information 
10¢0—1993 Constant equipment related data 
2000=Z232— Campaign identifying data 

S000 = So. Infrequent meaSurement data 

4000-4999 Space vehicle related data 

5000-5997 Weather data 

COlG—6Go 57 Real-time solution data 

TOUG= ao So Measurement data 

S000—8999 Post analysis data 


The data type identifiers of the Data Item List specify the 
orijer and the format of how the data is to be read. 

Data used in the analysis were read from the user 
exchange tapes using a slight modification of the sample 
input program contained in reference 19 (App. A). 
keal-time solution information (data type 6010) was read 
from the tapes and written to mass storage and retained. 
Other information such as measurement data types (data type 


920), real-time solution configuration (data type 1150), and 
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antenna location (data type 2210) were read manualiy fron 


the first page of the printouts of the tapes. 


C. COORDINATE TRANSFORMATIONS 


The coordinates of the stations used in the field 
testing were given in NAD-27 geodetic coordinates. The 
real-time soiution output of the GEOSTAK Receiver are in 
WGS-72 Cartesian coordinates. To make comparisons of the 
given positions to the observed positions, the given posi- 
tions were transformed to WGS-72 Cartesian coordinates. The 
transformation was done using the Abridged Molondensky 
formulas contained in the NGS‘'s conversion program (D035-41) 
for the HP-41CV programmable calculator. The progran 
prompts the user for the values of the origin shift between 
Coemeatatrer systems. ine x, Y, and 4 Origin shifts for 
Cratstonming [lem NAD=279to WGS-72 are -22, 157, and 176 
meters, respectively [Ref. 20]. These values differ fron 
the origin shift values supplied with the transformation 
program by NGS. The NGS values are -12.547, 155.804, and 
iWioc 109, for the X, Y, and Z origin shifts, respectively. 
The program generates Semimajor axis and inverse flattening 


parameters of both reference systems (Tab. II). 


TABLE II 
Transformation Program Parameters 


Bilipsoid a (meters) de 
Clarke 1866 6,378,206.4 294.9786982 
WGS 72 6,378,135.0 298.2600000 
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D. SAUVELIITE UPDASES 


The ephemeris elements and clock model updates occurred 
at identical times according to the data recorded in 
field tests. An approximation Ofetne Update title's ca cl. 
lated by subtracting the Age of Data (AODE) from tae ephen- 


eris reference time (t ods LeeCay 


update time = aa = 7 AQODE Cap 
The update times changed by -496, 1052, or 1552 seconds when 
the computation of equation 3.1 was performed tor each 
Satellite with subsequent occurrences of the broadcast 
ephemeris on the recording. The change is due to packing of 
bits in the transmissions of the AODE values. An approxiia- 
tion of the update times is made, within a few minutes, by 
using the first AODE values that change to smalier numbers 
(on the order of i to 2 hours) from previously recorded 


large values (on the order of 24 hours). [{Ref. 14] 


Ee AVERAGING METHODS 


Relative position determinations were made using reai- 
time solutions from two different stations at the same GPS 
time tags. Computer programs were written to read the real- 
time solutions from mass storage and compute the discrepan- 
cies between the observed and the given WGS-72 coordinates 
of pairs of stations at common GPS times. The beginning and 
ending times of data acguiSition varied amony the three 
receivers operated during the same time period; thus, the 
effect of comparing solutions at the same time tags resulted 
in the elimination of some of the real-time solutions fron 
the analysis. In effect, the tape with the earlier starting 
time waS spooled forward to correspond to the starting time 


of the first real-time solution recorded by the receiver 


ae 


that began operation at a sligntly later time. Similarly, 
the time of the last solution common to both receivers is 
the last solution recorded by the receiver that ended opera- 
tVOn MEL ESt 

A mean position was determined in two dirferent ways to 
Show the effect of the ephemeris and clock uvdates on the 
real-time solutions and to eliminate real-time solutions 
with large discrepancies from the averaging. In general, 
relatively large discrepancies are evident in roughly the 
Ticst hour of receiver operation. The real-time solution 
discrepancies of pairs of stations were averayed beginning 
with the values one kour after and two hours after the first 
solution time common to both receivers. In these tests, 
beginning the averaging one hour after the initial solutions 
corresponds to averaging the solutions that were recorded a 
short time before or after the update times and through the 
remainder of the recording session. Beginning the averaying 
two hours after the initial solutions corresponds to using 
only those values that were recorded much after the uvdates 
and through the remainder of the session. Real-time solu- 
tions are treated as independent observations in the statis- 


tical analysis. 


Fe. COMPUTATIONS 


The distances between stations ({d), in the WGS-72 rerer- 
ence system, are computed by equation 3.2 where xX, ,Y¥,, and 
4, ave the WGS-/2 Cartesian coordinates of one of the 


stations, and Xie y and Z., are the coordinates of the 


OY 
other station of the side. The side-lenygths are calculated 





a= Vix, - x)? + (Y¥,- ¥,)2 + (2, - 2.) (ann 


Zo 


for both the given coordinates and the mean values of the 
GEOSTAR real-time solution coordinates. The discrepancies 
between the two are computed by subtracting the given 
distances from the observed distances. Discrepancies are 
used in the averaging to save computer time and to reduce 
LOWnd-OfE Crror. 

The absolute value of the point position discrepancies 
are computed by eguation 3.3 where the subscripts o and g 


represent the observed and given values, respectively. 





5 eG. - 1.)2 + 1, = 39) eee (3.3) 


The standard error of a single observation was calcu- 
lated as a measure of the precision of the real-time solu- 
tions. The means and standard errors of a single, 12-second 
real-time solution were calculated for the coordinate 
discrepancies and the point position discrepancies. 
Comparisons are made for the two averaging methods (Sect. 
III E).~ The mean WGS-72 Cartesian coordinate discrepancies 
were used to compute the final observed positions. The 
observed positions are computed by subtracting the discrep- 


ancies from the given coordinates. 


G. GEODETIC COMPARISONS 


Geodetic accuracy standards for horizontai control 
reguire certain relative accuracies between directly 
connected adjacent points. The distances between points are 
computed along the eilipsoid - the geodesic. To compare the 
accuracies of the observed positions with the standards for 
geodetic control, the WGS-72 Cartesian coordinates of the 
observed positions are transformed to NAD-27 geodetic coor- 
dinates. The transformations were done using a modification 
of NGS's D035 program for the HP-41CV programmable 
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calculator. The modified version uses the Same parameters 
PseiOseeeceontatnecaninetne Orlginal program (fab. iI). 

The geodetic inverse computation is used to calculate 
the geodetic distance and azimuths between two known points. 
The inverses between given poSitions and between observed 
positions were computed and compared. 

To simulate the poSitionirg of an unknown statior, one 
station of a side is held fixed and the other station is 
located using the direct computation. The coordinates of 
the known position and the distance and azimuth to the 
unkrown position are the input variables and the coordinates 
of the unknown station are computed. The distance and 
azimuth calculated by the inverse computation of the 
observed point poSitions are used in this calculation. The 
direct and inverse computations were done on the IBM 3033 
based on Puissant's Coast and Geodetic Survey Formula 
[Ref. 21). 

The standard errors of the discrepancies of the geodetic 
position determinations for entire recording sessions are 
calculated aS an approximation of the repeatability. To do 
so, the two point poSition discrepancies of each station on 
days 61 and 68 must first be averaged to determine a single 


value for the station. 


2D 


IV. RESULTS ANDePESCUSSTON 


a = = = == = ee eee A SS SS 


A. GIVEN POSITIONS AND DISTANCES 


Data in this report are from relative position tests of 
Macch 1, 8, and 9, 1984. MTIhree receivers were operated 
Simultaneously over a period of approximately six hours. 

The antennas of the March Ist tests were located over a 
second-order position and two of its reference marks. The 
triangle has side-lengths of approximately 30 metcEes-. Qine 
tests of March 8th and 9th were conducted at first-order 
Survey points separated by distances of 14 to 25 kilometers. 
The data of the tests at one of the positions on March 9th 
waS not available at the time of this writing. The data 
gathered on March ist will be defined as the short base line 
tests and the 14 to 25 kilometer testing will be termed 
medium base line tests. Tables III and IV give the names 
and the NAD-27 coordinates oz the stations and Table YV gives 
the distances and azimuths in each triangle. Distances and 
aZimuths were calculated using the inverse computation 
(SCCt. (ebEr G)-. 

The NAD-27 coordinates of the given positions were 
transformed to WGS-72 Cartesian coordinates uSing a NGS 
transformation program (D035-41). The eliipsoidal heights 
are corrected for the heights of the antennas above the 
Station marks prior to the transformations. The neights of 
the antennas above the stations (data type 2110) were read 
manually from the first page of the printouts of the tapes. 
The distances between the points in WSS-72 Cartesian coordi- 
nates were calculated using equation 322 tal. Vie 

The computed side distances for the GS-72 Cartesian 


coordinate system are in all cases greater than or equal to 
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the Computed geodetic distances. In the medium base iine 
tests, the differences were found to increase by .0035 
TABLE iif 
Station Names 
J.D. Number Name Abbreviation 
61 85019 Oo oS. DTA oS 
Bo0 20 Poteet O5603 Reg) DMA RM 
B02 1 Nike Osos. Ray 2° OMA RM 2 
637569 go022 Beeson BGERO 
sa 027 BpLOt KNOE 1335 PaO TL 
*” 85028 BALD KNOB 1935 BNOQB 
TecertMOneNGte Weel 2n Narch 9th analysis 
TABLE IV 
NAD-27 Station Coordinates 
Station North Latitude West Longitude Ht. (m.) 
JMR S0SR 22 a6 "274255 O70 ees oe 9 7558 Zo 55 
RM 1 30022 57 54052 970M aS 367531041 Zo 
Ew 3085-22357". 06234 OT Otte 55.0 1009 235.26 
BZERO 50° - 73m 6 to 0238 Cg Oats 8205790 238.4 
PILO® B00 OSs Ot 3587 So Oe 28 238i ..03922 Daeg 
ENB 300 -7ere 0 7 01000 OF S35 5 Os) 700 204.4 
percent of the length of the line. Differences of .911, 
-748, and .534 meters were computed for the side-lergths of 


ApOEOXxEMately 25.5, 22.06, 


These ditfterences occur because the Cartesian distances are 


and 14.3 kilometers, respectively. 


computed between each point; the geodetic distances are 


computed from the projections of the points on the ellipsoid 


which does not account for the ellipsoidal heights. The 


distances between BZERO and PILOT are equal to the nearest 


TABLE V 


NAD-27 Geodetic Distances and Azimuths 


sae Distance (n.) AZimuth 
JMR - RM1 26.255 34 00) die uits 7 9570 
IMR - RM2 27059 699 06' 37"%.54898 
RM1 — Rf!2 37.759 1139 09' 19%.14446 
BZERO = PILOT 25 Sip cae 175° 32° welees537s 
RZERO - BNOB 14,346. 955 1139 51' 26". 80840 
BNOR - PILOT 22,556. 303 2099 40! 09",53915 
TABLE VI 


WGS-72 Side-length Distances 


J.D. side Distance (n.) 
61 JMR - BM1 Per. 255 
JMR - R42 27.593 
RM1 - R&2Z 3d. / 58 
68 BZEPO - PILOT Za, 495. | oo 
BZERO = BNGE 14,347.489 
BNOB =] 2sme7 Leeson se Oot 
69 BZERO - PILOT 25,493.188 


millimeter for tests run on days 68 and 69 even though the 
heights of the antennas above the stations were different on 
both days. 
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Be) OA Gee Le UPDATE TLNES 


The ephemeris and clock update tines occurred approxi- 
Mately ome hogr aiter the start of data acquisition in these 
mesSts.) “iene abe appboximately 25 Co 30 recordings of the 
broadcast ephemeris (data type 411)) during the six-hour 
recording periods. The upjiate times were computed (egqn. 
3-1) using the AODE values that chanyed from large numbers 
to much smaller numbers in subsequent recordings of the 
Lroadcast ephemeris (Tab. VII). Satellites with PRN codes 
9 and 11 were updated approximately one-half hour after the 
Mercessonr Coses Go ana won days ol™@and 69. Acsording “to 
the recorded iata, all four satellites were updated at 


approximately the same time on day 68. 


TABLE VII 
Satellite Update Times (Seconds) 








woe Ne COdes toe AODE to. - AODE 
6 1 6,8 S76,000 6, 144 S71, 056 
OG il Ser 0.0.0 4,096 S757 0 4 
E68 6,8,9,17 378,900 6,144 Selo 0 
69 ee 460,800 4,096 456,704 
oe 464,400 6,744 Oe. 90 


Cee oknnvrOr eon THE RERL-TINE SOLUIIONS 


The real-time solutions were read fron the GPS User 
Tapes and retained. The solutions of pairs of stations, 


with time tags of real-time solutions common to both, were 


Ze 


used in the analysis. The disckepameies Serre on tire 
observed coordinates and the yiven WGS-72 coordinates 
(observed - given) were computed and plotted over tine 
(Figs. 4.1 - 4.8). The vertagcal@anrows 1aeens EEaumes 
indicate the times of satellite updates (Tab. VII). The 
tick marks on the horizontal (time) axis represent one hour. 
In general, the X, Y, and 2 coordinate discrepancies 
begin with large values in the initial solutions and then 
become smaller and more consistent after satellite updates 
and over time. There is a Significant improvement in the 
discrepancies of the observed coordinates near the time of 
the update of PRN codes 6 and 8 at station JMR and near the 
update of codes 9 and 11 at station R'%1 on day 61 (Figs. 
G.1and 4.2). The precision of the disctCepancies in each 
coordinate increase at all three stations on day 61, a short 
time after the update of codes 9 and 11. The discrepancies 
at stations BZERO and PILOT on day 68 decrease approximately 


one hour after the initial solutions (Figs. 4.4 and 4.5). 


20 \GSai2) Point Postre 


Discrepancies 


— 


The absolute values or the point position discrepan- 
cies of pairs of stations with common time tags were 
computed (ecn. 3.3) for five-minute time segments and 
plotted with the standard errors over time (Figs. 4.9 - 
4.16). in general, the point position discrepancies and the 
standard errors decrease after satellite updates and over 
time. Except for the general trend, there is no obvious 
cocrelation of the accuracies (discrepancies) and the preci- 
Sion of the five-minute time segments at all stations on the 


Same day. 
3. WGS-72 Side-length Discrepancies 


Felatively high accuracy side-lenyth determinations 


acre evident at certain times; however, the times of these 
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STATION PILOT 
J.D. 68,1984. 
MEDIUM BASELINE TESTS 
5 MINUTE SEGMENTS 
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solutions are not consistent at all stations during the same 
Gay (Figs. 4.17 - 4.23). The discrepancies in the lengths 
of sides JME-RM2 and RM1-RM2 are relatively small near time 
tay 382632; however, the side-length discrepancy of JMR-RM1 
at approximately the same time is relatively larger than 
meer CetermindtLions ©f that side om day 61. Ihe side- 
length discrepancies of JMB-kM2 and RM1-RM2 are -. 044 and 
037 meters at time tags 382632 and 382600, respectively. 
mnie disecrevaney Of JNR-RM1 at time tag 382600 1s .955 
meters; whereas, the discrepancy is .476 meters at time tag 
386400. Therefore, the most accurate sije-lenjth determina- 
tions do not necessarily occur at approximately the sane 
time for all three sides. 

The most accurate side-length determinations do not 
necessarily occur when the magnitudes of the point position 
discrepancy vectors of two stations are nearly the sane 
because of differences in the directions of the discrepan- 
cies. The side-length discrepancies determined on day €5 
POGetne side BZERO-BNOB allustrate this point (Fig. 4.21). 
Th2 point position discrepancies for stations BZERO and BNOB 
_are 10.889 and 10.816 meters at tine tag 377808 and the 
side-length discrepancy is -.151 meters. At time tag 
384408, the side-lenyth discrepancy is -.036 meters and the 
magnitudes of the point position discrepancies are 9.295 and 
14.053 meters. The side-length discrepancy is smaller. at 
the later time tag even though the difference in the magni- 
tujes of the point position discrepancies are jreater at 
that time. 


D. AVERAGING RESULTS 
1. Determination of the Method of Averaging 


Two averaging schemes were used to arrive at mean 


position discrepancies (Sect. III E). The means and 
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Figure 4.18 Point Position and Side-length Discrepancies 
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Figure 4.19 Point Position and Side-length Discrepancies 
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Figure 4.22 Point Position and Side-length Discrepancies 
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Figure 4.23 Point Position and Side-length Discrepancies 
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standard errors of discrepancies were computed for tne X, Y, 
and Z coordinates and the absolute values of the totai point 
POET On GiSscCrEepancmescei(Tabs. Viiti@and If). 

The standard errors of the point position discrepan- 
cies for the solutions that were averaged beginning two 
hours after the initial solutions are, in general, smaller 
than those that were determined using ali soiutiorns begin- 
ning one hour after the initial solutions. The average 
standard error of the absolute value of the point position 
giscrepancies for all pairs of determinations is 2. 302 
meters using the former averaging method, compared to 1.422 
meters uSing the latter averaging method. Much of the 
cifference between these two standard errors can be attrib- 
uted to the relatively large point position discrepancies 
computed for stations BZERO and PILOT on day 69 between the 
ficst and second hours after the initial solutions (Figs. 
i> did | OpeeteglectingeSstatzons BZERO and PILOT on day 
69, the average of all Standard errors of the point position 
discrepancies are 1.€57 meters and 1.333 meters for the two 
different averaging methods. 

The elimination of the real-time solutions recorded 
between one hour and two hours after the start of data 
acjuisition results in smaller standard errors of the point 
position discrepancies in the averaginy process. The 
remainder of the analysis uses the mean coordinates deter- 
mined by the averaging method in which only the real-time 
solutions recorded two hours after th2 initial solutions are 
weed SiTa b.c 9 ee. 


z. WGS-72 Results 


The Y coordinate discrepancies are larger than the 
discrepancies of X and Z coordinates at all stations with 
the exception of those at station BZERO on day 69. The mean 


observed X coordinates are consistently differeat from the 


a0 
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given X coordinates by approximately 5 to 6 meters at all 
stations on all days of data acquisition having saneaverage 
standard error of .564 meters for all determinations. The Z 
coordinate has smaller discrepanemes at Aosta tons moneaa 4 
68 and 69 compared to the other coordinates. fhe average 
standard errors of the Y and Z coordinate discrepancies ior 
all determinations are 1.600 and 10560 meters mrcspec iver, 
(Tape). 

The meah point position discrepancies in the wWGS-72 
Cartesian coordinate system range between 16.650 meters at 
station RM2 on day 61 and 8.197 Meters at station 3S 7eo son 
day 69. The average of all the mean point position discrep- 
ancy determinations is 12.685 meters. The standard errors 
of the mean point position discrepancies range between .501 
meters at station BZERO on day 69 and 2.802 meters at 
station BNOB on day 68, with an average value of 1.383 


meters for all determinations (Tab. IX). 


E. GEODETIC COMPARISONS 


The observed WGS-72 cCoonpdinates Of Gach Sere romeo 
computed by subtracting the observed discrepancies from the 
given coordinates. The observed W3S-72 Cartesian coordi- 
hates are then transformed to NAD-27 geodetic coordinates 
using a modified version of NGS's 0035-41 transformation 


program. 


1. Discrepancies of Point Position Determinations 


The point poSition discrepancies of all pairs of 
determinations are computed for th2 NAD-27 coordinates (Tab. 
X). AzZimuths are measured from north. The average point 
position discrepancy is -.08336 seconds of latitude, =. 260388 
seconds of longitude, and +9.391 meters in ellipsoid height. 


The horizontal position discrepancies of the observed point 
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positions for all determinations range between 6.498 meters 
at station PILOT on day 69 and 8.530 metcES at Stacren eae 
on day 68. The average distance discrepancy along the 
ellipsoid is 7.596 meters. The azimuths from the given 
positions to the observed positions range between 129° 41° 
58".04469 at station BZERO on day 68 and 84° 48! 17%.17181 
-at station PILOT on day 69. The obServed positions have an 
average azinuth of 109° 25° S625630 fcomee ess ten 


positions. 
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2. Discre 


There are two determinations of point position 
discrepancies on days 61 and 68 as a result of comparing 
pairs of stations in the computations. an estimate of the 
repeatability of an entire recording session Can be calcu- 
lated by first averaging the two point position discrepan- 
cies for each station then computing the standard error of 
the eight determinations. 

The average discrepancies of entire recording 
sessions are -.07980 seconds of latitude and -.25903 seconds 
of longitude with standard errors ef 9039 7eeani e077 oc 
seconds, respectively. The average height discrepancy 1s 
8.658 meters above the station with a standard error of 
4.748 meters. The average horizontal distance jliscrepancy 


is 7.538 meters with a standard eCEbomprot 3714 meegers. 
3. Relative Position Comparisons 


The distances and azimuths between pairs of stations 
are computed for both the given coordinates and the observed 
coordinates by the inverse computation. The differences of 
tw are computed by subtracting the given distances and 
azimuths from the observed values (Tab. XI). Six of the 


seven observed side-lengths are shorter tnan the distances 


60 


TABLE ss 


Discrepancies of Observed Distances and Azimuths 


J.D. side Dist. (m.) Azimuth 
6 1 
JRM-RM1 1.397 -44" 44.2009 
JMR-RA2Z -0.006 “40 77* 4&O",. 3464 
RM 1-RM2 -0.231 -19 14* 30",.5102 
5 8 
BZERO-PILOT - 3.128 -1". 24259 
BZERO-BNOB -0.428 -26" .09039 
BNOB-PILOT -0.774 97 .49565 
59 
BZERO-PILOT -4.556 -1", 05927 


computed from the given coordinates. In the short base line 
tests, the differences in azimuth are of the order of 
degrees compared to differences of the order of seconds of 
arc in the medium base line tests. In general, the azinuth 
differences decrease as the side-length differences increase 


in both the short base line and nedium base line tests. 


4. kelative Position Discrepancies 


Relative positioning is used to determine an unknown 
position relative to the position of a known station. The 
dicect computation calculates the coordinates of a position 
Given the coordinates of a known position and the distance 
and azimuth to the unknown poSition. The distance and 
azimuth between pairs of observed positions are used in the 
Pee dee ederons wikeand IRM] are taken as fixed in the 
Short base line tests and the poSitions of RM1 and RM2 are 
calculated. Stations BZERO and BNOB are used to locate 


stations BNOB and PILOT in the medium base lin2 tests. The 


61 


discrepancies of the calculated positions are computed in 
terms of latitude, longitude, elevation, distiagecweazenae 
and relative accuracy between stations (Tab. XII). The 
elevation discrepancies are determined by computing the 
difference of the height discrepancies between pairs of 
point position solutions. 

Felative position determinations diiffer from the 
given positions by an average of 1.452 meters in the short 
base line tests and 2.712 meters in in the medium base line 
tests. The relative accuracies between stations range 
between 1:13 for the side JMR-RM71 on day 61 andi 1:17,418 for 
the side BNOB-PILOT on day 68. The average relative accu- 
racy of the medium base line tests is 1:9,712. Station 
PILOT was determined to relative accuracies of 1:8,139 and 
1:5,593 with respect to station BNOB on subseguent days. 

The data acquired in the short base line tests 
provide useful information in the absolute position anal- 
yses; however, this scenario 1s not a practical application 
of relative position determinations. It is not justifiable 
to measure distances of the order of thirty metars with such 
expensive equipment when the distances can be measured accu- 


rately with steel tape or other inexpensive equipment. 
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V. CONCLUSIONS AND RECOMMENDATIONS 

Based on the results of the medium base line tests, the 
real-time solutions of the TIG100/7SEOSTAR ave Caparo ot 
establishing relative positiops to Thim=d-order, (@8a-2 i: 
standards. Two receivers, one locatei at a known position 
and the other separated by 14 to 25 kilometers, can locate 
horizontal control in real-time. @ This assumes [ive Neuse: 
visibility of four satellites and a means for the outputeor 
a solution based on the averaging of real-time solutions 
recorded at 12-second intervals. 

It is found that relative positioning using GPS is more 
accurate than absolute point positioning. The average 
Gistance discrepancy of the absolute position determinations 
is 7.538 meters compared to relative position averades of 
1.452 meters in the short base line tests and 2.712 meters 
in the medium base line tests. Tha absolute point position 
determinations have a one Sigma repeatability of .714 meters 
in distance discrepancy. The relatively high repeatability 
of absoiute point position determinations resulted in rela- 
tive poSition accuracies of 1:9,712, on average, tor 
distances of 14 to 26 kilometers. 

Third-order, Class IE aceciitacies Gan be obtagmenuwi a 
the GEOSTAR Feceiver in real-time. A typical operational 
unit consists of two, 2-man teams; each team supplied with 
the GEOSTAR field set. One team sets an antenna over a 
known position and the other team locates at the point to be 
surveyed. Both receivers should gather data simultaneously 
during the periods of satellite visibility (during Phase iI 
testing). The teams should agree to begin the averaging 
algorithm at a predetermined time tag, preferably a short 
time after the uploads to the satellites. Satellite upload 


times can be requested beforehand from Vandenburg AFB. 
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Additional averaging methods should be explored. [In 
regard to the real-time solutions, averaging analyses could 
ke performed to determine the amount of data acguisition 
time required to achieve certain accuracies. An analysis of 
CHhiS Nature wwould bLewenmportart for Spesial military applica- 
tions where an accurate position 15 reyuired in less than 
five hours of receiver Operation. 

In the Fall 1984 the GEOSTAR will undergo hardware and 
SObUWonemNOdIrlcation oeroOre acull-scals production. It is 
planned that the software will be nodified to include real- 
time solutions of short time-interval observations; these 
solutions will be different from existing real-time solu- 
tions based on a cumulative least squares method. [{Ref. 22] 
Felative position determinations using the real-time solu- 
tions should improve with the new software because the solu- 
tions of two receivers can be compared simultaneously with 


no influence from preceding solutions. 
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GPS USER-TAPE READ-WRITE PROGRAM 


CK Fe a I ee oe A A Ee I a eK OK 


OueG) ©) [O): Gee) sO Gee) FO Or Ome) One) Oa sOe to 


THIS PROGRAM IS A SLIGHTLY MODIFIED VERSION OF THE 
EXAMPLE INPUT PROGRAM CONTAINED IN "A STANDARDIZED 
EXCHANGE FORMAT FOR NAVSTAR GPS GEODETIC DATA" BY 
BY V. DAN SCOTT AND J. GARRY PETERS, APPLIED RESEARCH 
LABORATORIES, THE UNIVERSITY OF TEXAS AT AUSTIN, MARCH 
10, 1983. THIS PROGRAM CAN BE EASILY CHANGED TO READ 
AND WRITE ANY DATA TYPE ON GPS USER-TAPES SUPPLIED BY 
NSWC. TO WRITE ALTERNATIVE DAIA TYPES, CONSULT THE 
ABOVE PUBLICATION TO DETERMINE THE DATA TYPES THAT 
CORRESPOND TO THE DATA OF INTEREST. EXAMPLE: TO PRINT 
OUT EPHEMERIS INFORMATION (DATA TYPE 4110), CHANGE 
THE GO TO STATEMENT FOLLOWING THE IF STATEMENT THAT 
TESTS FOR THIS DATA TYPE BELOW TO TRANSFER CONTEOL 
TO WRITE STATEMENTS WITH THE PROPER FORMAT. THE 
FORMATS OF ALL DATA TYPES ARE GIVEN IN THE ABOVE 
PUBLICATION. THE JCL FOR THIS PROGRAM IS FOR AN 
TBM 370 MODEL 3033. THE TAPES ARE RECORDED WITH 1600 
BPI, A LOGICAL RECORD LENGTH DF 80, A BLOCKSIZE OF 
6400, AND ASCII CODED. 


Co FE KK A ee FR a ee ke eo oo eo KK KK KK KK KKK KK K 


c 


//RAKOW JOB (1643,0130),*RAKOWSKY',CLASS=F 
//*MBRIN ORG=NPGYM1.1643P 


77 


EXEC FORTXKCLG 


Jf BOR Coe Ds ax 


c 


Cc 
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Mino eel DAG ge > LGT 


PVeroeeeele ly PRE, CONT R,NRED, ENDEX,DIYPE,JUNK,ASTER, 
clink nN i, Pius, CHARE,C HAKRF, CHART, NRID 
Rei a tk, O LG GZ, 5S KIP 


Ue Oy eye Uo” + /PCHARE/. &*/, CHARF/* FY /, 
tecunn/ hve phA wks */ ,SENT/'>°/ 


DIMENSION ICONTR(19) 


100 CONTINUE 
READ (1,10) ICN,IREP,ICONTR 
10 FOFPMAT(A1,13,19I4) 


120 CONTINUE 
IF (ICN.EQ. ASTER) GO TO 150 
IF (IGN-£O.PLUS) GO TO) 100 
TF (GCN. BO. CHARE) GO To 100 
IF (ICN.EQ.CHARF) GO To 500 
IF (ICN. EQ.CHARI) GO TO 500 
IF (ICN.29.SENT) GO T) 500 


150 CONTINUE 


TPE(ZEEP. EQ 20) ~LREP=1 
NREP=0 
ICN=BLANK 


200 Conn evs 


IN DEX=0 

NREP=NREP+1 

TP (DREPSLIaG)) GO TO 300 

IF (NREP.GT.IREP) GO TO 100 
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309 


GO TO 300 


CONTINUE 


IF (ICN.NE. BLANK) GO TO 120 


INDEX=INDEX+1 


IF (INDEX.GT.19) GO TO 200 
DTYPE=ICONTR (INDEX) 
IF (DTYPE.EQ.0) GO TO 200 


TP (DTYPE. 80.910) GO Tomo 
IF(DIYPE.EQ.920) GO TO 201 


IF (DTYPS. EQ. 1010) 
IF (DTYPE. EQ. 1020) 
IF (DTYPE.EQ. 1030) 
IF (DTYPE. EQ. 1040) 
IF (DTYPE.EQ. 1050) 
IF (DTYPE. EQ. 1110) 
IF (DTYPE. EQ. 1150) 
IF (LTYPE. EQ. 2010) 
IF (DTYPE.EQ. 2020) 
IF (DTYPE. FQ. 2030) 
if (DIYPE. BOs 200} 
IF (DTYPE. EQ. 2120) 
IF (DTYPE.EQ. 2130) 
IF (DTYPE. EQ. 2140) 
IF (DTYPE. EQ. 2210) 
IF (DTYPE. EQ. 2250) 
IF (DTYPE.EQ. 2270) 
IF (DTYPE. EQ. 3010) 
IF (DTYPE. EQ. 3210) 
IF (DTYPE. EQ. 3220) 
IF (DTYPE. EQ. 3250) 
IF (DTYPE.EQ.40 10) 
IF (DTYPE. FQ. 4050) 


GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 


-GO 


GO 
GO 
GO 
GO 
GO 
GO 
GO 
GO 


TO 
TO 
TO 
0 
TO 
TO 
TO 
oe 
TO 
Le 
ie 
TO 
EO 
TO 
TO 
a 
TO 
TO 
To 
EO 
ALLO, 
LO 


m 
4 


204 
205 
204 
205 
204 
204 
202 
20 1 
201 
20 1 
203 
20 1 
201 
20 1 
202 
20 1 
20 1 
202 
20 1 
20 1 
201 
201 
201 
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IF(DTYPE.EQ.4110) GO TO 204 
IF (DTYPE.EQ.4120) GO TC 202 
LF (DTYPE.EQ.4210) GO TO 203 
TEI E. EOs 1 OS0\e GO TO 202 
IF (DTYPE.EQ.4270) GO TO 202 
Pa DEYPES EO. S010) GO TO 201 
IF (DTYPE.EQ.6010) GO TO 6010 
IF (DTYPE.EQ.6020) GO TO 202 
IF(DTYP&.E9.6100) GO TO 201 
IF (DTYPE.EQ.6110) GO TO 203 
PET LPE. PO. 70 a0) sce TO 207 
IF (DTYPE.EQ.7110) GO TO 202 
IF(DTYPE.EQ.7120) GO TO 201 
IF (DTYPE.£Q.7210) GO TO 201 
IF (DTYPE.£0.7220) GO TO 201 
Bas cO. 7250) Go TO 201 
TED T YES. Bem 2 LOCO NTO™ 202 
IF (DTYPE.£Q.7250) GO TO 202 
IF (DTYPE.E0.7260) GO TO 203 
WRITE (4,31) 
31 FORMAT('DATATYPE ERROR") 


201 CONTINUE 
READ (1,11) ICN 
11. FORMAT(A1) 
IF ((ICN. NE. BLANK) .AND. (ICN.NE.PLUS)) GO I) 120 
SOLO, 300 


202 CONTINUE 
READ(1,11) ICN 
IF ((ICN.NE.BLANK) .AND. (ICN.NE.PLUJS)) GO IO 120 
READ (1,11) NCN 
GG TO 390 


Zoo. COMI ins 
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READ (1,11) ICN 

IF ((ICN.NE.BLANK) AND. (ICN.NE.PLUS)) GO TO 120 
READ(1,11) NCN 

READ(1,11) NCN 

GO TO 300 


204 CONTINUE 
READ (1,11) ICN 
IF ((ICN.NE.BLANK) .AND. (ICN.NE.PLUS)) GO TO 120 
READ(1,11) NCN 
READ(1,11) NCN 
READ(1,11) NCN 
GO TO 300 


205 CONTINUE 
READ(1,11) ICN 
IF ((ICN. NE. BLANK) .AND. (ICN.NE.PLUS)) GO Td 120 
READ(1,11) NCN 
READ(1,11) NCN 
KEAD(1,11) NCN 
READ (1,11) NCN 
GO TO 300 


60 10 ,GCiT Inve 


FEAD(1,12) ICN ,NRID, Boner ey 

IF ((ICN. NE. BLANK) .AND. (ICN.NE.PLIS)) GO ITD 120 
READ (1,13) (NEN, 2, SiG? ,ShGmaolGlsiG7> skne 
WRITE (4, 14) DD Ne kee 

WRITE (4,15) SGX SiG oes 


C READ FORMATS 


12 FORMAT(A1,12,8X,D22.15,D22. 15,D17. 10) 
13. FORMAT (A1,2D17.10,4E10. 3} 


ie 


©, -§G212 ORMATS 


18 “S@RMAT (D22.15, P22. 15,D17. 10,017.10) 
Sere A (31%) e10 20x, £1023 >4% 6 10. 3) 


c 
GO TO 300 
20 PCONTENUE 
SOF. 
END 
Cc 
C 


//30-FTOIFO01 DD UNIT=3400-3,VOL=SERK=mytape, 

VEE ISP=(OLOP PASS LAB Sk=t2,NL,, LN), 

// DCB=(RECFH=FB, LR ECL=80, BLKSIZE=6400, DEN=3, DP TCD=Q) 
V7iaOetnoenoO)l DD SDSN=NSslS1643.FLLENAME, UNIT=3330V, 


// DISP= (NEW, CATLG,DELETE), MSVGP=PUB4B, 
// UCB= (RECFM=FB, LRECL=80, BLKSIZE=649)) 
/* 

as 
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APPENDIX 38 


PROGRAM TO COMPUTE REAL-TIME SOLUTION DISCREPANCIES 


Ck a ok kkk kk ak kkk kkk kkk kok ac a ak ak ok ke oka oko dk ak kok ok ea ok 


©) Gi -@) GO) OO) BGG). OO) 2) eG) Oe) Oe ee) 


THIS PROGRAM COMPUTES THE DL FFERENCE SOP ee oa 
AND Z COORDINATES, THE TOTAL POINT POSTITTON eter 
ANCIES, AND THE SIDE-LENGIA DISCREPANCIES 2) haa 
STATIONS AT COMMON GPS TIMES. THiS PROGR ee ee 
WRITES THE ABOVE DATA BEGENN ING Wilde i ae ieee oe 
TIME SOLUTIONS OF EACH TAPE, THEN SPOOlS siete. 
THE EARLIER RECORDING SESSION [CO COnnsSseonp Oe. 
GPS TINE JF THE TAPE THAT STARTED Con Oi NGS eae ee 
TIME. A CONVERSION OF THE GIVEN GaODETTC Cl 7.7 aes 
TO WGS-72 CAKTESIAN COORDINATES MUST Pikstes. Delle. 
THESE VALUES ALONG WITH COMPUTED SiDE Ueto t rere Ce 
ARE ENTERED BELOW FOLLOWENG THE TYPE DEP iit ioe 2 ae 
OUTPUT IS IN THE FORM: “GIVEN VALUES = REAL eee 
SOLUTION Valine 


CK KK HE EK eK i eo ee te eo ae te eK ee ee ok ie ee oko ok ae ek eK Ke ok kk 


Cc 


//RAKOW JOB (1643,0130),"RAKOWSKY',CLASS=B 
//*MAIN ORG=NPGVMN1.1643P 


1/ 


EXEC FORTXCLG 


J/EORTS Syst N Dies 


Cc 


REAL*8 TTAG1, TTAG2,X%1,X2, 71,712,721, o2, DELX, Deine 
*¥ DIST,DVALUE, SIGK1,SIGXK2, SiG vig siGy 275i Ga ome ae 
* DELDST,X1IGIVN,X2GIVN,YIGIVN, Y2GIVN,Z1GIVN, Z2GIVN, 
* XIDIF, X2DIF, YIDEF, Y2DIF, Zio, 220nr, OFS ef pip ereseuee 


DVALUE=25493.188 


2 


30 


X IGIVN=- 740577.548 
Y 1IGIVN=-5456773.662 
Z1GIVN=32076609 .081 
X2GIVN=-740336.501 
WeZGrVN==596970 72020 
Z2GIVN=3185698.432 


READ (1,10, END=500) ITAG1,X1,¥1,21 
Raa (a iiypesiGx?,SiGY1,SiGZ 1 
HELD (2 ho, END= 500) TRAG2,X%2,1 2, 22 
PED 25 19) SIG Kl, SiGY?, SiGZ2 


TEriiet en To TTAGZ)) COM TOM 


CONTINUE 

Pee Oj gmiennG2.k25.02 -7 2 
POR oe MEST G © Shey? ,SIGZ2 
IF (TTAG2.GE.TTAG1) GO TO 30 
GO TO 2 


CONTINUE 

READ(1,10) ITAG1,X1,¥1,21 
READ; lEorGMIeSTGyY 1,Sicz1 
Die 2AGweGherrAce) GO TO 30 
GO TO 3 


GONTINUVUE 

PERT PAG aGiswLAGZ), GO TO 2 

PC RAGZeGlettAGI) GO TO 3 
DELX=X1-X2 

DELY =f ta e2 

DELZ=Z1-2Z2 

Pistols so DEL XarZ2+DELY**2Z+DELZ** 2 } 
DEUS CiSslt—=DVALUE 


i 


OQ? 


©) 


X1DIF=X1i-X1GIVNW 
Y WDLEPSY I~ Y 16a 
ZVOT YSZ 1=Z VGN 
X2DIF=X2-X2GIVN 
Y2DiIb=Y2—¥ 2Gna 
Z2DI F=Z2-Z2GIVN 
OFSETI=DSQORT (X IDL E**2 tye ee 
OFSETZ=DSQRT (X2DIF**2+ Y2DIF** 2+Z2D IF ¥¥ 2) 


WRITE (4,42) TTAG1,TTAG2 

WRITE (4,43) X1DIF,Y1DIF,Z1DIF 
WRITE(4,41) OFSET1 

WRITE (4,43) X2DIF,Y2DIF,Z2DIF 
WRITE(4,41) OFSET2 
WRITE(4,41) DELDST 


READ (1,10, END=500) GEaGaeele Ihe 
PEAD (1,11) SIGX1,SIGY1)sueom 

READ (2,10, END=500) TTAGOH Ko er 2 
READ (2,11) SiGX2, SIG’ secure 

GO To 30 


READ FORMATS 


10 
is 


PORMAT (022.15, D22. 15, Dilealeeson 7 eet on 
FORMAT (31X,D17. 10, 9x, UG ore 3) 


VE EEE SeeORMATS 


41 
42 
4 3 


500 


FORMAT (D22. 15) 
FORMAT (2D22.15) 
FORMAT (3D22. 15) 


CONTINUE 
S51 OF 
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END 


e 
//30.FTO1FOO01 DD DSN=MSS.S1643.MYTAPE, 
We DISP= (OLD, KEEP) , MSVGP=PUBU4S 
//30.FTO2F001 DD DSN=MSS.S1643.FILENAME, 
// DISP= (OLD, KEEP) ,MSVGP=PUB4B 
COnEtdsn001 DDEDSN="Ss5.5 1643-CHHCK, UNIT=3350, 
/* DISP= (lizW, KEEP) ,VOL=SER=MVS004,SPACE=(CYL, (4,4)), 
/% DC B= (RECFM=FB, LRECL=80, BLKSIZE=6400) 
/* 
//30.FTOUFOCT DD DSN=MSS.S1643. filename, UNIT=3330V, 
// PEoP=(NEN GAT LG, DELETE) ,MSVGP=PUBUB, 
Ve DC B= (RECFM=FB, LRECL=80, BLKSIZE=6400) 
a 
VG 
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AODC 
AODE 
Polgtsss ets 
Aoerien 
CDU 

DMA 
DMAHTC 
GEOSTAR 
GUT 

GPS 
Mbps 
NGS 
NAD-27 
NAVSTAR 
NSWC 
NSHC/DL 
PPS 

PRN 

SPS 


WGS-72 


APPENDIX C 

LIST OF ABBREVEATIONS 
Age of Dita (Clock 
Aje of Data (Ephemeris) 
Applied Research Laboratories: Univ. of Texas 
American Standard Code for Enio. Pnterenange 
Control and rsSpoa sore 
Defense Mapping Agency 
Defense Mapping Agency Hydro./Topo. Center 
Geodetic Satellite Tracking and Ranging 
Greenwich Mean Time 
Global Positioning Systen 
Million bits per second 
GPS Master Control Station 
North American Datum of 1927 
Navigation Satellite Tracking and Ranging 
Naval Surface Weapons Center 
Naval Surface Weapons Center/Dahlgren Labs. 
Precise Positioning Service 
Pseudo-Random Noise 
Standard Positioning Service 


World Geodetic System of 1972 
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